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Abstract 

We present preliminary results of a maximum-likelihood Dalitz-plot analysis of the charmless 
hadronic decays to the final states ^it^tt^ and K^tt^tt^ from data corresponding to an inte- 
grated on-resonance luminosity of 166 fb _1 recorded by the BABAR experiment at the SLAC PEP- 
II asymmetric-energy B Factory. We measure the total branching fractions B{B^ — > tt^-k^tt^ 1 ) = 
(16.2±2.1±1.3)xl0~ 6 and B{B ± K^n^ir^) = (61.4±2.4±4.5) x 10" 6 , and provide fit fractions 
and phases for intermediate resonance states. 
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1 Introduction 



The study of charmless hadronic B decays can make important contributions to the understanding 
of CP violation in the Standard Model as well as to models of hadronic decays. The B-meson 
decay to the three-body final state can proceed via intermediate resonances formed from two of the 
particles. These two-body states can interfere with each other and with the nonresonant three-body 
decay. The three-body state is unique in the search for weak phases as it is possible to determine the 
strong phase variation for overlapping resonances. Studies of these decays can also help to clarify 
the nature of the resonances involved, not all of which are well understood. A full Dalitz-plot 
analysis is necessary to correctly model this interference and extract branching fractions. 

Observations of B-meson decays to the ir^ir^ir^ and K^tt^tt^ three-body final states have 
already been reported by the Belle and BABAR collaborations using a method that treats each 
intermediate decay incoherently [1, 2, 3]. Belle has also reported a preliminary Dalitz-plot anal- 
ysis of the decay B ± — ► K ± 7r T ir ± [4]. We present a preliminary Dalitz-plot analysis for both the 
B ± — > 7r ± 7r =F 7r ± and B ± — > K ± it T it ± decay modes. 

2 The BABAR Detector and Data Sample 

Here we present preliminary results from a full amplitude analysis based on a 166 fb _1 data sam- 
ple containing 182 million BB pairs collected with the BABAR detector [5] at the SLAC PEP-II 
asymmetric-energy e + e~ storage ring [6] operating at the T(4S) resonance at a center-of-mass en- 
ergy of y/s = 10.58 GeV. An additional total integrated luminosity of 16 fb _1 was recorded at an 
energy 40MeV below this energy and is used to study backgrounds from continuum production. 
The charm decay B + — ► D°tt + , D° — ► K + ir~ is used as a calibration channel as it has a relatively 
high branching fraction. 

Details of the BABAR detector are described elsewhere [5]. The specific components used for 
this paper are charged particle tracking provided by a combination of a silicon vertex tracker 
(SVT), which consists of five layers of double-sided detectors, and a 40-layer central drift chamber 
(DCH) in a 1.5-T solenoidal magnetic field. This allows a transverse momentum resolution for the 
combined tracking system of (J Pt /pt = 0.0013j>t © 0.0045, where the sum is in quadrature and pt 
is measured in GeV/ c. Charged-particle identification is provided by combining information on the 
average energy loss (dE/dx) in the two tracking devices and the angle of emission of Cherenkov 
radiation in an internally reflecting ring-imaging Cherenkov detector (DIRC) covering the central 
region. The dE/dx resolution from the drift chamber is typically about 7.5% for pions. The 
Cherenkov angle resolution of the DIRC is measured to be 2.4 mrad, which provides nearly 3a 
separation between charged kaons and pions at a momentum of 3 GeV/c. 

3 Event Selection and Reconstruction 

B-meson candidates are reconstructed from events that have four or more charged tracks. Each 
track is required to have at least 12 hits in the DCH, a minimum transverse momentum of 100 
MeV/c, and a distance of closest approach to the primary vertex of less than 1.5 cm in the trans- 
verse plane and less than 10 cm along the beam axis. Charged tracks identified as leptons are 
rejected. The B-meson candidates are formed from three-charged-track combinations and parti- 
cle identification criteria are applied. The average selection efficiency for kaons in our final state 
that have passed the tracking requirements is ~ 80% including geometrical acceptance, while the 
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misidentification probability of pions as kaons is below 5% at all momenta. The kaon veto on pions 
in our final state is ~ 98% efficient. The I?-meson candidates' energies and momenta are required 
to satisfy appropriate kinematic constraints, as detailed in Section 5. 

4 Background Suppression and Characterisation 

The dominant source of background comes from light quark and charm continuum production. 
This background is suppressed by imposing requirements on event-shape variables calculated in 
the T(4S) rest frame. The first discriminating variable is cos 9t, the cosine of the angle between 
the thrust axis of the selected B candidate and the thrust axis of the rest of the event. For 
continuum background the distribution of | cos 9t\ is strongly peaked towards unity whereas the 
distribution is uniform for signal events. We require |cos6*t| < 0.575 for B + — > 7r + 7r~7r + and |cos6*r| 
< 0.65 for B+ -► K+7t~7t + . 1 

Additionally, we make requirements on a Fisher discriminant T [7] formed using a linear com- 
bination of nine variables representing the angular distribution of the energy flow of the rest of the 
event into each of nine two-sided concentric 10° cones around the thrust axis of the reconstructed 
B [8]. 

Other backgrounds arise from BB events. There are four main sources: combinatorial back- 
ground from three unrelated tracks; three- and four-body B — > DX decays, where X represents 
other particles in the final state; charmless four-body decays with a missing particle and three-body 
decays with one or more particles misidentified. In the case of charm decays with large branching 
fractions these backgrounds are greatly reduced by vetoing the appropriate region of the two-body 
invariant-mass spectra. The rejected decays and the invariant-mass veto ranges are given in Table 1. 

Table 1: The invariant-mass veto ranges (in GeV/c 2 ) for intermediate resonances for 7r + 7r~7r + and 
K + ir~ir + . The leptons in the J/tp and i[)(2S) decays are misidentified as pions. 



Resonance 


7T + 7T 7T + 


K + ir ir + 


J/0 -► t+f 

ip{2S) -► e+e- 

D° -► K+TT- (or 7T+7T-) 


3.05 < m n+7V - < 3.14 
3.64 < m n +„- < 3.73 
1.70 < m I+r < 1.93 


2.97 < m w+w - < 3.17 
3.56 < m n + n - < 3.76 
1.80 < m K+7T - < 1.90 



The remaining charm backgrounds that escape the vetoes and backgrounds from charmless B 
decays are studied using a large sample of Monte Carlo (MC) simulated BB decays equivalent to 
approximately five times the integrated luminosity for the data. Any events that pass the selection 
criteria are further studied using exclusive MC samples to estimate reconstruction efficiency and 
yields. The ttt-es an d Dalitz distributions of the BB backgrounds, which are used in the likelihood 
fits, are then normalised to the total number of predicted BB events in the final data sample. For 
■k + it~it + , we expect 200 ± 14 5-related background events, dominated by the decays B + — > Kgir + 
and B + — ► K + it~ir + . For K + it~it + , we expect 315 ± 17 background events and the dominant 
backgrounds come from J>-meson decays to states containing a D or D* and a p or 7r, the decays 
B° — > p + K~ and B + — ► r/K + , and the nonresonant decay B° — ► K + 7t^tt°. 

1 Charge- conjugate states are implied throughout this section. 
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A further background in this analysis comes from signal events that have been misreconstructed 
by switching one or more particles from the decay of the signal B meson with particles from the 
other B meson in the event. The amount of this background is estimated from MC studies, and for 
both B + — ► 7r + 7r~7r + and B + — > K + ir~ir + is found to be a very small effect that accounts for less 
than 2% of the final data sample in the signal box (defined in Section 5). As such it is neglected 
in the analysis. 

Both the continuum and S-related backgrounds are modeled in the Dalitz amplitude fit using 
linearly interpolated 2-dimensional histograms. 

5 Final Data Selection 

Two kinematic variables are used to select a final data sample. The first variable is AE = 
E* B — \/s/2, the difference between the center of mass (CM) energy of the l?-meson candidate 
and \fs/2, where y/s is the total CM energy. The second is the energy-substituted mass tttes 

= (s/2 + pi • pB) 2 /Ef — p 2 B where p# is the B momentum and (Ej,pj) is the four-momentum of 
the initial state. The mean of the AE distribution is shifted by — 8.3MeV as measured from the 
calibration channel B + — ► D ir + . For tt^tt^-k^ we require —68.3 < AE < 51.7 MeV; for i^ ± 7r T 7r ± 
the requirement is —38.3 < AE < 51.7 MeV where the lower edge is tightened by 30 MeV to reduce 
K^tt^tt^ specific BB backgrounds. 

We define three regions in the AE-ities plane, illustrated in Figure 1. The signal box is defined 
by 5.271 < ttjes < 5.287 GeV/c 2 and events in this region are used in the amplitude analysis. 
The signal strip is defined by 5.20 < tties < 5.29GeV/c 2 and is used to determine the fraction of 
signal and qq continuum events in the signal box. A sideband area below the signal box, defined 
by 5.20 < ?rtEs < 5.26GeV/c 2 , is used to obtain the distribution of the qq continuum events in the 
Dalitz plane. 

We accept one I?-meson candidate per event in the signal strip. Fewer than 3% of events have 
multiple candidates and in those events one candidate is randomly accepted to avoid bias. 

After all selection criteria are applied the average efficiency for reconstruction of phase space 
7r ± 7r =F vr ± and i^ ± 7r =F vr ± MC events in the signal box is (13.00 ± 0.04)% and (13.26 ± 0.03)%, 
respectively, where the errors are statistical only. The efficiency across the Dalitz plot is uniform 
except for very small decreases near the boundaries. In the amplitude analysis, this is taken into 
account by calculating the efficiency as a function of position in the Dalitz plot using phase space 
MC. 

The mES distribution for events in the signal strip is used to determine the number of qq and 
signal events in the signal box. The signal component is modeled by a double Gaussian function 
with parameters obtained from phase-space MC. These parameters are fixed except for the mean of 
the core Gaussian. The qq continuum background is modeled using the experimentally motivated 
ARGUS function [9] with the endpoint fixed to the beam energy while the shape parameter is al- 
lowed to float. Finally, the BB background is modeled with an ARGUS function plus a Gaussian to 
account for peaking BB backgrounds. All parameters of the BB component, including the amount 
of peaking and nonpeaking BB background, are obtained and fixed from the MC simulation. The 
fraction of signal and qq events is allowed to float. Figure 2 shows the m,ES projections of fits to 
the data for both -k^-k^-k^ and K ± 7r T 7r ± . The \ 2 per degree of freedom for these projections is 
0.83 (1.13) for 7r ± 7r =F 7r ± (K ± n T n ± ). For 7r ± 7r =F 7r ± the total number of events in the signal box is 
2407; for K ± tt^tt ± the total number of events in the signal box is 3174. The extracted fractions of 
signal, qq continuum and BB backgrounds are given in Table 2. 
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Figure 1: AE-tues plane, showing the signal strip (entire populated region), sideband and signal 
box regions, as defined in the text, for — ► K^it^it^ 1 . Similar regions apply for B^ — > tt^-k^^. 
The events shown are B ± -» K± ■k^t: candidates from the data sample. 




Figure 2: toes distribution, together with the fitted PDFs: the data are the black points, the lower 
solid red area is the qq component, the middle solid green area is the BB background contribution, 
while the upper blue line shows the total fit result. All errors shown are statistical only. The left 
hand plot is for tt^tt^tt^ 1 ; the right hand plot is for K^ir^ir^. The vertical dashed blue lines show 
the signal region extent, while the dotted black line indicates the upper edge of the sideband. 



6 Dalitz Amplitude Analysis 

In terms of a Dalitz-plot analysis of the £>-meson decay to the final state hhh (where h = it or K) 
a number of intermediate states contribute and the total rate can be represented in the form: 



dT 



dm^dm^ 



l 13' 



™23) 



(1) 



11 



Table 2: Total number of events in the signal box, for -k^tt^tt^ and K^tt^tt^ candidates. Also 
shown are the fractions for each hypothesis calculated from the fit to the mES distribution for 
events in the signal strip. The errors arise from the combination of the fit errors on the signal 
and qq components, and the error on the BB background due to the uncertainties on the various 
branching fractions and MC efficiencies. 





7r ± 7r =F 7r ± 




Total Events 


2407 


3174 




Fraction of Events 


Signal 


0.159 ± 0.021 


0.467 ± 0.021 


qq background 


0.758 ± 0.018 


0.430 ± 0.014 


BB background 


0.083 ± 0.010 


0.103 ± 0.006 



where m 2 3 and m| 3 are the invariant masses squared of pairs of final state particles. In the case 
of ir^ir^-K^, m\ 3 = (p n ± +p 7rT ) 2 and = (p n ± +p^) 2 ; for i^" ± 7r =F 7r ± , m\ 3 = (p K ± +p 7rT ) 2 and 
m 23 = (Pir± +Pttt) 2 - The amplitude for a given decay mode is Cie l6i Fi{m\ 3 , «i| 3 ), where c$ and 
9i are the unknown real parameters of each partial decay mode, while Fi describes the dynamics 
of the amplitudes. These Fj consist of a product of the invariant mass and angular distribution 
probabilities: 

Fi = Ri(m) x Ti{cos9 H ) (2) 

where Ri{m) is the resonance mass distribution and Tj (cos Oh) is the angular probability distribu- 
tion. The angle Oh is defined as the angle between the momentum vector of one of the resonance 
daughters in the resonance rest-frame and the momentum vector of the resonance in the B rest- 
frame. 

To fit the data in the signal box, we define an unbinned likelihood function for one event to 
have the form shown in Eq. 3. The fit is performed allowing the amplitude magnitudes (cj) and 
the phases (0j) to vary. 



r(rn 2 2 , n f f x I YaU Cie i0 *Fi(ml 3 , mj 3 )| 2 e(m 2 3 , m 2 3 ) 



where 



+ fBB 



I IE l =iQe^F(m 2 3 ,m| 3 )| 2 e(m 2 3 ,m| 3 ) dm\ 3 dm\ 3 
Q(m 2 13 ,ml 3 ) 



B(rn\ 3 ,m\ 3 ) 



• m\ 3 and m| 3 are the invariant mass-squared values of the daughter pairs (1,3) and (2,3); 

• iV is the number of resonant and nonresonant contributions to the plot; 

• Fi is the dynamical part of the amplitude of the resonant or nonresonant contribution i; 

• Cj and 0i are the real parameters to be determined (— ir < 0i < n); 
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e(m1 3 ,m,2 3 ) is the reconstruction efficiency defined for all points in the Dalitz plot; 
Q(mf 3,77133) is the distribution of qq continuum background; 



• B(mf 3 , T7t2 3 ) is the distribution of BB background; and 

• fqq and f B g are the fractions of qq continuum and BB background events, respectively. They 
are determined from the ttjes fit and MC, respectively, and are fixed in the amplitude fit. 

The first term on the right-hand-side in Eq. 3 corresponds to the signal probability density 
function (PDF) multiplied by the signal fraction (1 — f qq - — f BB )- Since we can always apply a 
common factor to both the numerator and demoninator of the signal PDF, this analysis will only 
be sensitive to relative phases and magnitudes, and hence it is possible to fix the magnitude and 
phase of one component (p°(770) for 7r ± 7r =F 7r ± and IT" (892) for K ± 7t t 7t ± ). 

As the choice of normalisation, phase convention and amplitude formalism may not always be 
the same for different experiments, fit fractions are presented instead of amplitude magnitudes to 
allow a more meaningful comparison of results. The fit fraction is defined as the integral of a single 
decay amplitude squared divided by the coherent matrix element squared for the complete Dalitz 
plot as shown in Eq. 4. 

Fit Fraction = T f^ ^777^ 9' 9 , 9 , 9 • 4) 

J IEiCie^Fi(mf3,m^3)| 2 dm| 3 dm^3 

Note that the sum of these fit fractions is not necessarily unity due to the potential presence of 
net constructive or destructive interference. 



7 Physics Results 

7.1 B ± — ► T^Tr^Ti-i Results 

The nominal fit is performed with the resonances p°(770), p°(1450), /o(980), /2(1270), and a 
uniform nonresonant (NR) contribution. The masses and widths of the resonances are fixed to 
their world average values [10]. In this fit the p°(770) is fixed to have a magnitude of 1 and its 
phase is set to 0, since this is the dominant contribution to the Dalitz plot, and this choice reduces 
the statistical uncertainties for the other fitted components. We model all the resonances using 
relativistic Breit-Wigner lineshapes with Blatt-Weisskopf barrier factors [11] except for the /o(980), 
which is modeled with a Flatte lineshape [12] (to account for its coupled-channel behaviour due 
to the fact that it can decay to 7r + 7r - or K + K~). The nonresonant component is assumed to be 
uniform in phase space. The preliminary results of the nominal fit to the 2407 events in the signal 
box can be seen in Table 3, along with the total branching fraction (BF) and the average efficiency 
across the Dalitz plot weighted by the fitted signal distribution. 

Figure 3 shows the mass projection plots for the nominal fit. The four resonant contributions 
plus the single uniform phase-space nonresonant model are able to adequately describe the data 
within the statistical uncertainties. 

Further fits are performed to the data by removing one two-body component at a time from the 
nominal model; the results are shown in Table 4. Removing the p°(770), p°(1450) or ^(1270) com- 
ponents give significantly worse fit results. The omission of the /o(980) and nonresonant amplitudes 
give values of the remaining fitted components that are close to their nominal values. 
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Table 3: Results of the nominal fit for the — * tt^-k^tt^ mode. The first errors are statistical, 
while the second errors are systematic and are described in Section 8. 



Average Efficiency (%) 


13.0 ±0.9 ±0.6 


Total BF (xlO _ti ) 


16.2 ± 2.1 ± 1.3 


Component 


Fit Result 


p u (770) Fit Fraction (%) 
p°(770) Phase 


58.2 ±2.9 ±6.0 
0.0 (Fixed) 


p°(1450) Fit Fraction (%) 
p°(1450) Phase 


13.6 ±2.8 ±2.0 
±0.59 ±0.39 ±0.17 


/ (980) Fit Fraction (%) 
/o(980) Phase 


2.0 ±1.3 ±2.8 
±2.45 ±0.61 ±0.19 


/ 2 (1270) Fit Fraction (%) 
/ 2 (1270) Phase 


14.3 ±2.0 ±1.8 
-2.69 ±0.33 ±0.17 


NR Fit Fraction (%) 
NR Phase 


4.2 ± 2.0 ± 1.4 
±0.61 ±0.56 ±0.19 




Figure 3: Projection plots of the fit result for B — > n n^ir onto the mass variables m^^- and 
m^-i^i. The data are the black points with error bars, the lower solid red histogram is the qq 
component, the middle solid green histogram is the BB background contribution, while the upper 
blue histogram shows the total fit result. All errors shown are statistical only. The large dips in 
the spectra correspond to the vetoes in Table 1. 



Recent experimental results from e + e collisions at BES [13] show evidence of a low-mass tt + tt 
pole in data for J/ip — > ujir + 7r~ , known as the a. Analysis of data from the E791 experiment for 
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D + — > 7r + 7r~7r + [16] show similar results. Also a large concentration of events in the 1 = S-wave 
7T7T channel has been seen in the region around 500 — 600 MeV in pp collisions [14]. This pole is 
predicted from models based on chiral perturbation theory [15], in which the resonance parameters 
are M - iT/2 = [(470 ± 30) - i(295 ± 20)] MeV. 

Consequently the a resonance is predicted in S ± — ► 7r + 7r~7r ± decays. For this Dalitz-plot 
analysis the a resonance is modeled using the parameterisation suggested by Bugg [17]. The result 
of the Dalitz-plot fit is shown in Table 5. It can be seen that the inclusion of the a slightly changes 
the contributions of the other resonant components and also affects the phase of the nonresonant 
amplitude, but all results are consistent with those from the nominal fit, within the statistical 
uncertainties. 

Table 5 also shows the preliminary results of the fit to the data when we add the /o(1370) and 
XcO resonances to the model. Both contributions are found to be negligible. 

Table 4: Results of the fit to the Dalitz plot for — ► tt^tt^tt^ candidates, for the nominal fit 
and for fits performed with a different component omitted in turn from the nominal fit. All errors 
shown are statistical only. 





Nominal 


No p u (770) 


No p°(1450) 


No / (980) 


No / 2 (1270) 


No NR 


(— ln£) — (— hxC (nominal)) 




94.3 


8.0 


1.7 


18.9 


2.8 


p u (770) Fit Fraction (%) 


58.2 ±2.9 




71.2 ±3.0 


58.4 ±2.6 


66.4 ±3.8 


58.9 ±2.5 


p u (1450) Fit Fraction (%) 
p°(1450) Phase 


13.6 ±2.8 
±0.59 ± 0.39 


46.2 ± 3.5 
±0.59 (Fixed) 




14.8 ±3.0 
±0.58 ± 0.39 


16.1 ±3.6 
±0.30 ±0.41 


16.9 ±2.9 
±0.46 ±0.36 


/ (980) Fit Fraction (%) 
/ (980) Phase 


2.0 ± 1.3 
±2.45 ±0.61 


21.9 ±2.9 
±2.58 ±0.29 


4.5 ± 1.3 
±2.92 ±0.39 




1.6 ± 1.8 
±2.27 ±0.83 


1.5 ± 1.0 
±2.78 ±0.59 


/ 2 (1270) Fit Fraction (%) 
/ 2 (1270) Phase 


14.3 ±2.0 
-2.69 ±0.33 


23.2 ±2.6 
-2.95 ±0.24 


16.7 ± 1.4 
-2.09 ±0.22 


14.1 ±2.1 
-2.75 ±0.35 




14.8 ±2.0 
-2.84 ±0.33 


NR Fit Fraction (%) 
NR Phase 


4.2 ±2.0 
±0.61 ± 0.56 


34.2 ±3.7 
±1.84 ±0.23 


11.8 ± 2.5 
±0.82 ±0.39 


3.1 ± 1.7 
±0.27 ±0.58 


6.2 ± 2.7 
±0.07 ± 0.56 
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Table 5: Results of the fit to the Dalitz plot for B — > tt tt^tt candidates, for the nominal fit and 
for fits performed with a different component added in turn to the nominal fit. All errors shown 
are statistical only. 





Nominal 


With cr 


With / (1370) 


With xco 


(— hxC) — (— ln£ (nominal)) 




-3.9 


-0.1 


0.0 


/9 U (770) Fit Fraction (%) 


58.2 ±2.9 


52.8 ±6.4 


58.4 ±3.7 


58.1 ± 5.2 


p u (1450) Fit Fraction (%) 


13.6 ±2.8 


10.6 ±3.2 


13.6 ±2.8 


13.6 ±3.0 


p°(1450) Phase 


±0.59 ±0.39 


±0.88 ± 0.46 


±0.58 ±0.38 


±0.60 ±0.39 


/ (980) Fit Fraction (%) 


2.0 ± 1.3 


5.9 ± 1.6 


1.7 ±1.8 


1.9 ± 1.3 


/ (980) Phase 


±2.45 ±0.61 


±2.19 ±0.37 


±2.32 ± 0.79 


±2.45 ±0.61 


/ 2 (1270) Fit Fraction (%) 


14.3 ± 2.0 


12.4 ±2.5 


14.2 ±2.1 


14.2 ±2.3 


/ 2 (1270) Phase 


-2.69 ±0.33 


-2.44 ± 0.36 


-2.69 ±0.33 


-2.69 ±0.33 


NR Fit Fraction (%) 


4.2 ±2.0 


1.0 ±2.6 


3.8 ±2.7 


4.1 ±2.0 


NR Phase 


±0.61 ±0.56 


±2.99 ± 1.06 


±0.69 ±0.71 


±0.61 ±0.56 


a Fit Fraction (%) 




25.7 ±9.1 






cr Phase 




-1.62 ±0.21 






/o(1370) Fit Fraction (%) 






0.0 ±2.1 




/ (1370) Phase 






-1.7 ±2.2 




XcO Fit Fraction (%) 








0.2 ±7.5 


Xco Phase 








-0.7 ±4.7 
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7.2 B ± -> ^±71-^71^ Results 

The expected significant contributions to the K^-k^-k^ Dalitz plot can be identified from previous 
studies of this final state. Our model includes the following resonances: if* (892), if,* (1430), 
p°(770), /o(980), XcO and a nonresonant amplitude. The masses and widths of the resonances are 
fixed to their world average values [10]. In the fits we use if* (892) as the reference component 
and hence its magnitude and phase are fixed to 1 and 0, respectively. We use the same lineshapes 
described in Sec. 7.1 to model the resonance dynamics, except for if q (1430), where we use the 
LASS amplitude model. The dynamics of the Kir S-wave are not very well established and form an 
area of some disagreement within the community. Some favor the existence of the n pole [17], whilst 
others strongly oppose it. All, however, agree that there is strong evidence of resonant behavior 
around 1430 MeV, the ifo°(1430). The LASS experiment studied Kir scattering and as part of 
this study produced a description of the S-wave that consists of a resonant part, the if,* (1430), 
and an effective-range term [18, 19]. This amplitude is only measured up to around 2 GeV/c 2 in 
Kit mass, and so we curtail the effective-range term at the lower edge of the D° veto. Since the 
LASS amplitude contains both a resonant and nonresonant part, the results for if,* (1430) are not 
purely due to this resonance, but to the Kir S-wave as a whole. This model, with five two-body 
components plus a uniform nonresonant component, with the LASS amplitude for the if,* (1430) 
will be referred to as the "nominal" model. 

The nominal fit shows very good agreement with the data; a comparison can be seen in Figure 4. 
The preliminary results of the nominal fit to the 3174 events in the signal box can be seen in 
Table 6 along with the total branching fraction (BF) and the average efficiency across the Dalitz 
plot weighted by the fitted signal distribution. We have also used a relativistic Breit-Wigner or a 
Flatte lineshape for the if,* (1430) resonance with and without the addition of a k resonance, as 
suggested in Ref [17], to model the Kir S-wave but found the fits to be poor representations of the 
data. 

Table 6: Results of the nominal fit for the — ► K^ir^ir^ mode. The first errors are statistical, 
while the second errors are systematic, and are described in Section 8. 



Average Efficiency (%) 


12.8 ±0.8 ±0.5 


Total BF (xl0~ B ) 


61.4 ±2.4 ±4.5 


Component 


Fit Result 


if* (892) Fit Fraction (%) 
if* (892) Phase 


11.4 ±2.0 ±1.5 
0.0 (Fixed) 


if *°(1430) Fit Fraction (%) 
if *°(1430) Phase 


52.6 ±2.3 ±4.0 
±2.92 ±0.11 ±0.10 


p°(770) Fit Fraction (%) 
p°(770) Phase 


8.5 ± 1.9 ± 1.1 
±0.85 ±0.38 ±0.35 


/o(980) Fit Fraction (%) 
/o(980) Phase 


15.0 ±2.4 ±1.3 
-0.55 ±0.32 ±0.41 


XcO Fit Fraction (%) 
XcO Phase 


1.45 ±0.27 ±0.23 
±0.15 ±0.33 ±0.22 


NR Fit Fraction (%) 
NR Phase 


7.9 ±0.9 ±2.3 
+0.50 ±0.24 ±0.24 
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Figure 4: Mass pair projection plots for the nominal — ► if ^tt^tt^ fit. The left plot shows the 
K^it^ mass spectrum and the right plot shows the ir^ir^ mass spectrum. The data are the black 
points with error bars, the lower solid red histogram is the qq component, the middle solid green 
histogram is the BB background contribution, while the upper blue histogram shows the total fit 
result. All errors shown are statistical only. The large dips in the spectra correspond to the vetoes 
in Table 1. 



We refit the data removing each component in turn and find, in each case, that the fit worsens 
considerably. In some cases the parameters of the other components can vary dramatically. Table 7 
details the results of these omission tests. We also test for the possibility that there are further 
resonances in the Dalitz plot by repeating the fit using the nominal model but adding an additional 
component. The possible additional states we test for are the /2(1270), /o(1370) and /O (1450) 
resonances in the tttt mass spectrum and the iT| (1430), K (1680) and k resonances in the Kir 
mass spectrum. The results of these fits are shown in Table 8. In all cases the — In £ is slightly 
better than in the original fit. These checks are performed only on the B + sample. 



18 



Table 7: Results of the fit to the Dalitz plot for B + — > K + ir ir + candidates, for the nominal fit and for fits performed with a different 
component omitted in turn from the nominal fit. All errors shown are statistical only. 





Nominal 


No K* u (892) 


No if * u (1430) 


No p u (770) 


No / (980) 


No XcO 


No NR 


(— ln£) — (— In £ (nominal)) 




55.4 


156.9 


25.0 


79.9 


17.0 


21.5 


K* {) (892) Fit Fraction (%) 


9.5 ± 1.8 




19.9 ± 1.1 


9.8 ± 1.8 


10.3 + 4.0 


9.6 ± 1.7 


9.2 ± 1.5 


K^(U30) Fit Fraction (%) 
if *°(1430) Phase 


54.3 ±3.1 
+2.98 ±0.16 


64.5 ±9.6 
+3.03 + 0.34 




59.4 ± 3.8 
+2.96+0.16 


58.8 ± 3.9 
+3.03 + 0.15 


54.8 ±3.1 
+3.03 + 0.16 


54.2 ± 2.9 
+2.89 + 0.15 


p u (770) Fit Fraction (%) 
p°(770) Phase 


8.6 ± 1.9 
+ 1.48 ±0.45 


8.6+1.1 
+ 1.73 + 0.28 


13.1 ± 1.5 
+1.68 + 0.30 




14.4 ±5.0 
+ 1.55 + 0.57 


8.2 + 1.7 
+1.67 + 0.43 


11.0 + 2.5 
+0.45 ± 0.45 


/ (980) Fit Fraction (%) 
/ (980) Phase 


13.4 ±2.4 
-0.18 + 0.40 


13.3 + 2.0 
-0.18 (Fixed) 


32.3 ± 2.3 
+0.18 + 0.21 


16.8 + 3.0 
+0.06+0.40 




12.7 + 2.2 
+0.05 ± 0.39 


20.8 + 3.2 
-1.45 ±0.37 


XcO Fit Fraction (%) 
Xco Phase 


1.57 ±0.41 
+0.75 ± 0.47 


1.57 + 0.41 
+ 1.06 + 0.46 


1.46 ± 0.24 
-2.45 ±0.34 


1.57 + 0.41 
+0.67+0.47 


1.38 + 0.39 
+0.84 ± 0.48 




2.27 ± 0.56 
+ 1.08 + 0.47 


NR Fit Fraction (%) 
NR Phase 


9.5 + 1.3 
+0.71 ± 0.32 


7.35 ± 0.97 
+1.16 + 0.27 


43.2 ± 2.6 
-1.62 ±0.17 


13.7 + 1.8 
+0.66+0.32 


14.1 ± 1.8 
+0.99 ± 0.32 


12.2+1.5 
+0.89 + 0.30 





Table 8: Results of the fit to the Dalitz plot for B + — > K + it tt + candidates, for the nominal fit and for fits performed with a different 
component added in turn to the nominal fit. All errors shown are statistical only. 





Nominal 


With / 2 (1270) 


With / (1370) 


With /9 U (1450) 


With if* u (1680) 


With K 2 *°(1430) 


With k 


(— \n£) — (— In C (nominal)) 


— 


-2.1 


-2.1 


-1.9 


-0.7 


-3.7 


— 3.7 


if*" (892) Fit Fraction (%) 


9.5 ± 1.8 


9.0 ± 2.4 


9.5 ±2.2 


9.2 ±2.7 


9.2 ±1.9 


9.8 ±2.1 


9.4 ±2.1 


K£ u (1430) Fit Fraction (%) 


54.3 ±3.1 


53.5 ± 3.7 


53.7 ± 3.6 


55.5 ±4.6 


53.3 ±3.4 


50.7 ±3.1 


54.4 ± 3.6 


i^°(1430) Phase 


+2.98 ±0.16 


±3.02 ±0.16 


±3.00 ±0.16 


±2.92 ±0.17 


±2.97 ±0.17 


±3.00 ±0.16 


2.97 ± 0.16 


p°(770) Fit Fraction (%) 


8.6 ± 1.9 


7.7 ± 1.8 


7.7 ±1.7 


9.0 ±1.3 


8.5 ±2.0 


8.0 ±1.7 


7.1 ±1.6 


p°(770) Phase 


±1.48 ±0.45 


±1.47 ±0.45 


±1.66 ±0.44 


±1.24 ± 0.55 


±1.50 ±0.46 


±1.87 ±0.43 


±1.20 ±0.44 


/o(980) Fit Fraction (%) 


13.4 ±2.4 


13.8 ±3.5 


11.0 ±2.5 


14.3 ±4.1 


13.6 ±2.6 


13.3 ±2.7 


12.8 ±2.6 


/ (980) Phase 


-0.18 ±0.40 


-0.14 ±0.48 


±0.22 ±0.44 


±0.18 ±0.51 


-0.20 ±0.42 


±0.15 ±0.44 


-0.62 ±0.44 


Xco Fit Fraction (%) 


1.57 ±0.41 


1.53 ±0.41 


1.57 ±0.42 


1.61 ±0.46 


1.58 ±0.43 


1.55 ±0.43 


1.61 ± 0.42 


Xco Phase 


±0.75 ± 0.47 


±0.77 ±0.47 


±0.77 ± 0.47 


±0.54 ±0.49 


±0.80 ± 0.48 


±0.87 ±0.48 


±0.72 ± 0.47 


NR Fit Fraction (%) 


9.5 ± 1.3 


9.5 ± 1.3 


11.5 ± 1.6 


11.0 ±2.0 


9.0 ±1.3 


8.6 ± 1.5 


9.9 ±1.3 


NR Phase 


±0.71 ±0.32 


±0.75 ±0.31 


±0.74 ±0.32 


±0.53 ±0.38 


±0.73 ± 0.34 


±1.07 ±0.38 


±0.74 ±0.31 


/ 2 (1270) Fit Fraction (%) 




2.1 ± 1.9 












/ 2 (1270) Phase 




-0.09 ±0.49 












/ (1370) Fit Fraction (%) 






1.15 ±0.74 










/ (1370) Phase 






±0.38 ± 0.64 










p u (1450) Fit Fraction (%) 








6.5 ±2.5 








p°(1450) Phase 








-0.92 ±0.56 








K*°(1680) Fit Fraction (%) 










0.71 ± 1.24 






K*°(1680) Phase 










±2.35 ±0.85 






if 2 * u (1430) Fit Fraction (%) 












3.95 ±0.54 




if 2 *°(1430) Phase 












-0.35 ± 0.25 




k Fit Fraction (%) 














3.03 ± 0.89 


k Phase 














±2.90 ± 0.42 



8 Systematic Studies 



The charged particle tracking and particle identification uncertainties are 2.4% and 3.0%, respec- 
tively. There are also global systematic errors in the efficiencies due to the criteria applied to the 
event-shape variables (1.0%) and to AE and tues (1.0%). We also take into account the statistical 
uncertainty on the efficiency (due to the weighting by the amplitude model) for the individual and 
total branching fraction results (7% for B for — > K^ir^ir^). The uncertainty 

in the number of BB events is evaluated to be 1.1%. 

The systematic error on the efficiency variation across the Dalitz plot is calculated by performing 
a series of fits to the data where we vary the contents of each bin in the efficiency histogram according 
to binomial errors. This introduces an absolute uncertainty of 0.02 to 0.09 for the phases, and a 
fractional uncertainty of 0.4% to 4.9% for the fit fractions. However, for the average efficiency, and 
hence for the total branching fraction, this is a very small effect, evaluated as 0.1%. 

The systematic uncertainty introduced by the BB background and qq background has two 
components, each of which can potentially affect the fitted magnitudes and phases differently. The 
first component arises from the uncertainty in the overall normalisation of these backgrounds, whilst 
the second component arises from the uncertainty on the shapes of the background distributions 
in the Dalitz plot. The uncertainties on the fit fractions and phases due to the normalisation 
uncertainty are estimated by varying the measured background fractions in the signal box by their 
statistical errors. The maximum associated absolute uncertainty for the phase is 0.85 (0.05) for 
n^ir 11 ^ (K^ir^ir^) due to the qq background normalisation uncertainty and 0.03 due to the 
BB background normalisation uncertainty. These uncertainties are added in quadrature. The fit 
fractions are affected in a less uniform manner, with relative uncertainties in the range 0.7% to 15.4% 
(0.2% to 4.1%) for tt^tt^tt^ (K^tt^tt^). The uncertainties on the fit fractions and phases due to 
the Dalitz-plot background distribution uncertainty is estimated in the same way as the efficiency 
variation, namely varying the contents of the histogram bins in accordance with their Poisson 
errors. To be conservative, each phase for -k^-k^-k^ has been given an associated uncertainty of 0.1 
due to the qq background distribution uncertainty and 0.1 due to the BB background distribution 
uncertainty, which are then added in quadrature. For K^tt^tt^ the variations between modes are 
larger and each mode is therefore treated individually. The range of uncertainties in the phases 
is 0.06 to 0.36 due to qq background and 0.04 to 0.18 due to BB background, which are again 
added in quadrature. The fit fractions are affected in a less uniform manner, with the relative 
uncertainties ranging from 1% to 8% (1% to 25%) for tt^tt^tt^ (K^tt^tt^). 

Possible biases due to the fitting procedure are investigated using extensive Monte Carlo sim- 
ulations. We assign an absolute systematic error up to 3% for the fit fractions and 0.1 for the 
phases. 

9 Summary 

Tables 3 and 6 show the preliminary results from the nominal fits together with their statistical 
and systematic errors. 

For the decay B^ — > n^ir^n^ the nominal fit to the Dalitz plot is performed with the resonances 
p°(770), p°(1450), /o(980), /2(1270) and a uniform nonresonant contribution. The total branching 
fraction is consistent with the previously measured value [2, 10]. We can estimate the branching 
fraction of B^ — ► / o°(770)7r =l: by multiplying its fit fraction by the total branching fraction for 
B ± — > . We find this to be (9.4 ± 1.3 ± 1.0) x 10 6 , which is consistent with the previously 
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measured value [20, 10]. The removal of any of the fit components gives worse likelihood values, 
especially for the p°(770), p°(1450) and /2(1270) resonances. It is found that the /o(1370) and XcO 
resonances have no contribution to this Dalitz plot, while there is some evidence for a contribution 
from the a. 

For the decay — > K^ir^ir^ the nominal fit to the Dalitz plot is performed with the resonances 
K*°(892), i-C *°(1430), p°(770), / (980), XcO, a uniform nonresonant component and the LASS 
parameterization of the scalar component of the Kir spectrum. The total branching fraction result 
is consistent with the previously measured value [2, 10]. We can estimate the branching fraction 
of B ± -> K*°(892)ir ± ,K*°(892) -> K^tt? by multiplying its fit fraction by the total branching 
fraction for B ± —> K ± tt t tt ± . This yields a value of (7.0±1.3±0.9) x 10~ 6 , which is smaller than that 
reported by earlier analyses that do not fit over the full Dalitz region but is consistent with the Dalitz 
analysis reported by Belle [4]. The resonance behavior around 1400 MeV is successfully modeled 
with the LASS parameterization of the Kir S-wave. The removal of any of the fit components 
results in a significant worsening of the fit likelihood and a large change in one or more of the 
remaining amplitudes. Addition of further resonances (/ 2 (1270), / (1370), p°(1450), K|°(1430), 
K*°(1680) or k) does not cause a significant change in the fit likelihood, and in most cases the 
extra component is measured to be compatible with having zero fit fraction. In addition, their 
presence in the fit does not significantly alter any of the other parameters and consequently has no 
significant effect on the fit fractions of the original components. The fit is stable with respect to 
different parameterizations of the Flatte line shape of the /o(980) . 

Future iterations of the analysis will use the information from the separate B + and B~ sam- 
ples to make measurements of the charge asymmetries of the various intermediate decay modes, 
potentially allowing observation of direct CP violation. 
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